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ABSTRACT 

Several authors have shown that precise measurements of transit time variations of 
exoplanets can be sensitive to other planetary bodies, such as exo-moons. In addition, 
the transit timing variations of the exoplanets closest to their host stars can provide 
tests of tidal dissipation theory. These studies, however, have not considered the effect 
of the host star. There is a large body of observational evidence that eclipse times of 
binary stars can vary dramatically due to variations in the quadrupole moment of the 
stars driven by stellar activity. In this paper we investigate and estimate the likely 
impact such variations have on the transit times of exoplanets. We find in several cases 
that such variations should be detectable. In particular, the estimated period changes 
for WASP-18b are of the same order as those expected for tidal dissipation, even for 
relatively low values of the tidal dissipation parameter. The transit time variations 
caused by the Applegate mechanism are also of the correct magnitude and occur on 
timescales such that they may be confused with variations caused by light-time travel 
effects due to the presence of a Jupiter-like second planet. Finally, we suggest that 
transiting exoplanet systems may provide a clean route (compared to binaries) to 
constraining the type of dynamo operating in the host star. 
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1 INTRODUCTION 

Since the discovery of the first ex oplanet around a solar- 
like star bv lMavor fc Queloj ^99^), the exoplanet field has 
bloomed with over 400 currently known. While the major- 
ity of exoplanets have been discovered through radial ve- 
locity measurements of the Doppler wobble eff'ect, it is the 
systems that exhibit transits that are most highly prized. 
These planets are crucial for determining exoplanet bulk 
densities (via radii and mass measurements) as well as their 
atmospheric properties (via infrared measurements of their 
day /night variations and transmission spectroscopy). 

Transiting exoplanets, however, also ofi'er the opportu- 
nity to detect other planets within the system since an addi- 
tional planet may alter the period of the observed transits. 
This can occur in two ways. In the first case the gravitational 
infiuence of the perturbing body can can alter the orbital pe- 
riod of the transiting exoplanet directly. This effect is par- 
ticularly strong for planets in mean motion resonances and 
can even allow Earth-massed objects to be detected, while 
'exo-moons' orbiting the transiting planet itself also induce 
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a similar effect (e.g. ISimon et al.ll200'i1 ). In the second case, 
a perturbing mass in a wider orbit can cause the transiting 
planet / star system to wobble around the barycentre, again 
altering the observed transit times but by changing the light 
travel-time. 

Although searching for transit timing variations (here- 
after, TTVs) can p otentially uncover th e existence of Earth- 
mass objec ts (see [ Gibson et al.l Rabus et al 

] I2OO9I : 

iBcan 200^: iGibsoiTe t al. 2009; Miller-Ricci ct al. 20(3, and 
references therein for recent observational studies), there 
are other effects that can lead to TTVs. These include 
the p recession of orbits due to general relativistic ef- 
fects l|Pal fc KocsisI |2008| ). tidal dissipation, torques due 
to the spin-induced q uadrupole moment of the host star 
(|Miralda-Escudel |2002| ) . perturbations of transit times due 
to star spots, as well as reorientation of the planetary or- 
bit with respe ct to the Earth as a result of proper motion 
l|Rafikovl l2009l). However, there is a wealth of observations 
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I2OO8I ) that have shown quasi-periodic variations in eclipse 
times over timescales of years to decades that are compara- 
ble to, or larger than, the effects being searched for amongst 
transiting exoplanets. 

The favoured explanation for these observed variations 
in the orbital peri ods of eclip sing binary stars is known as the 
Applegate effect (jApplegatei 1992). This mechanism invokes 
magnetic activity cycles in the low-mass components of such 
binaries to redistribute angular momentum within the inte- 
rior of the star, thereby changing the stellar quadrupole mo- 
ment which lead s to changes in the orbital period o f the com- 
ponents. Later, iLanza. Rodono fc Rosnen ()l998l l proposed 
that the Applegate mechanism could also be driven by ef- 
fectively converting rotational kinetic energy and magnetic 
energy back and forth. Regardless of the details of the ex- 
act physical mechanism at work, the Applegate effect should 
also operate in most exoplanet systems since the host stars 
are (by selection) low-mass stars with a convective outer 
layer which should exhibit some form of dynamo activity. 
It will therefore be important to know the magnitude of 
the Applegate effect for exoplanet systems when interpret- 
ing any TTVs. In this paper we briefly review the Applegate 
m echanism in the n ext section, before applying the analysis 
of lApplegat3(ll992l ) to estimate the effects on known transit- 
ing exoplanet systems. Finally, we look at the implications 
that the Applegate effect has for TTV work in detecting ad- 
ditional planets as well as for the measurement of the level 
of tidal dissipation in very-hot Jupiter's. 



2 PERIOD CHANGES IN BINARY STARS - 
THE APPLEGATE EFFECT 

Many types of binary stars show evidence for changes in 
their orbital periods revealed most easily through eclipse 
times. If a star suddenly increases its orbital period P by 
an amount AP, then the eclipses will arrive progressively 
later and later until, after a time T, they are delayed by an 
amount At = TAP/P with respect to an ephemeris based 
upon the initial period. These variations can be tracked 
through comparison of observed times to those calculated 
assuming linear ephemerides via the so-called 'O ~ C dia- 
grams'. 

The period changes observed in close binary stars have 
typical magnitudes of AP/P ~ 10~^ (which can take on ei- 
ther sign), and thus significant deviations can build up over 
time. One famous example is the 2.87 day-period binary Al- 
gol, which has exhibited deviations from linearity of the or- 
der of 3 hours over the 200 years it has been followed. Cata- 
clysmic variable stars, detached white dwarf/main-sequence 
binaries, RS CVn stars and W UMa stars all exhibit similar 
variations on timescales of years to decades. These varia- 
tions, which are at best quasi-periodic, are more often than 
not too large to be explained by long-term effects such as nu- 
clear evolution of the stars, mass loss through winds or angu- 
lar momentum loss from gravitational radiation or magnetic 
braking. The latter in particular struggles when confronted 
wit h orbital periods that in cr ease a s well as decrease. 

lApplegate fc PattersonI |l983) realised that period 
changes in binary systems without either mass or angu- 
lar momentum loss could be driven by variations in the 
quadrupole moment of one or both stars. Taking just one 



star to have a quadrupole moment Q (in the equato- 
rial planej) and mass M then, as shown in equation 4 of 
Applegate (1992), its companion orbits in a gravitational 
potential of the form, 
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where r is the distance from the star. The orbital speed is 
given by v'^ = rdcj>/dr and therefore, from equation [T] if Q 
increases (the star becomes more oblate) then the gravita- 
tional field in the equatorial plane of the star also increases. 
In order to balance for the increased gravity, the companion 
then requires to increase its centrifugal acceleration n^/r at 
constant angular momentum {rv is constant). Thus v must 
increase and r must decrease and hence the orbital period 
decrease s. The opposite i s true if Q decreases. Under this 
scenario, lApplegatd (|l992l ) showed that the resulting period 
changes are given by (their equation 7): 
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where R is the radius of the star and a is the orbital sepa- 
ration of the components. 

The common denominator in all the binaries that ex- 
hibit period variations of the type described above is that at 
least one of the components is a low-mass star with a convec- 
tive envelope capable of sustaining a magnetic field generat- 
ing dynamo. Given that the variation timescales are of the 
order of years to decades, stellar activity cycles are the prime 
candidate for effecting the cha nges in the stellar quadrupole 
moment . Models of the 1980's (|Applegate fc Pattersonlll987l : 
I Warned [l988) supposed that dynamo-generated magnetic 
fields distorted t he star from its equilibr ium shape. These 
were criticised bv lMarsh fc Pringld l| 19901 ) since such distor- 
tions leave the pressure and gravity gradients unbalanced 
(the star is driven from hydrostatic equilibrium). This un- 
bala nced force needs t o be balanced by the magnetic field, 
and iMarsh fc Pringld (|l990l ) showed that the stars had in- 
sufficient luminosity to drive such changes on the observed 
timescales. 

This issue was quickly resolved by lApplegatd (|l992l ) 
where magnetic activity was still invoked to drive the Q 
variations but rather than forcing the star from hydrostatic 
equilibrium, the magnetic fields are supposed to drive an- 
gular momentum transfer within the star. For instance, if 
angular momentum is transported from the core to the en- 
velope of the star, the star will become more oblate overall. 
This required much less energy than the earlier models since 
the star performs a transition from one st ate of hydrostatic 
equilibrium to another. With this model lApplegatd (Il992l ) 
was able to explain the observed period changes in binary 
stars, and the Applegate effect has continued t o survive the 
test o f time apart from a few refinements (e.g. iLanza et al.l 
1 19981 ). 



3 ESTIMATING THE PERIOD CHANGES IN 
EXOPLANET SYSTEMS 

The same process that occurs in binary stars will also oc- 
cur in exoplanet systems in the case where the host star is 
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magnetically active. Since the majority of exoplanet hosting 
stars known to date are lower main-sequence stars with con- 
vective outer layers, these stars should harbour some form 
of magnetic field generating stellar dynamo. 

For any giv en energy budget (AE), the equations of 
lApplegatd (|l992l ) allow the magnitude of the period change, 



AP 

P ' 



(3) 



to be calculated numer ically which , for th e sake of complete- 
ness, we outline here. lApplegatd (|l992l ') adopted a simple 
stellar model in which a thin outer shell of mass Ms is spun 
up by addition of angular momentum AJ from the interior 
of the star. The energy required to do this is given by his 
equation 28, 



AE = fldrAJ + 
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21, 



(4) 



where Q.dr = Os — il* is the angular velocity of differential 
rotation between the outer shell {Q.s) and the stellar interior 
leff is the effective moment of inertia given by, 



(5) 



where I3 and 7* are the moments of inertia of the outer shell 
and the stellar interior, respectively. Typically, for an outer 
shell mass of Ms = O.IM, Is — I, and therefore in equation[4] 
we can substitute 27e// = Is, where Is — 2/3MsR^- 

The Qdr term on the left-hand side of equation |4] is 
normally small (but see later) and can be set to zero. We 
can therefore re-express equation [4] as, 

3(Ajf 

In the model outlined in lApplegatd |l992'), variations 
in the quadrupole moment, Q, of the star are driven by the 
stellar activity cycle. Magnetic fields are supposed to drive 
the angular momentum transfer from within the star. For 
instance, if angular momentum is transported from the core 
of the star to its envelope, the star will be come more oblate 
and its quadrupole moment will increase. lApplegatd (| 19921 ') 
computes the rate of change of the stellar quadrupole as a 
function of the angular momentum transport in his equation 
26, 
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where Q is the stellar angular rotation velocity. 

As outlined in Section [21 a change in the stellar 
quadrupole moment leads to a corresponding change in the 
orbital period given by equation[2l where a is the orbital sep- 
aration between, in this case, the star and the planet. We 
now have three equations, with equation [6] relating AE to 
A J, equation [7] relating AQ to A J, and finally equation [2] 
relating AP to AQ. Given that stellar activity cycles are 
rather variable, there is no well-defined amplitude or period 
of variation to be expected from Applegate's models. The 



one constraint that we can invoke is a restriction on the en- 
ergy budget allowed to drive the quadrupole moment. With 
this in mind, we can use the last 3 equations to relate the 
(as yet undefined) energy budget to the change in the orbital 
period AP giving. 
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If we assume that the power available to drive the 
quadrupole changes is some fraction, /, of the stellar lumi- 
nosity then we obtain a total energy budget of AE = fLT, 
where T is the timescale over which the quadrupole changes 
occur. Since stellar activ ity cycles tend to be of th e order 
of years or decades (e.g. ISaar fc Brandenb^I^ll999l ). there 
is potentially a large energy budget available for driving or- 
bital period changes. If this energy is supplied by the nuclear 
luminosity of th e star with no en ergy storage in the convec- 
tion zone then (jApplegatdll992h the star will exhibit RMS 
luminosity variations of 



ALr 



AE 



(10) 



We set the available energy budget such that the luminosity 
variations never exceed some fraction a of the total stellar 
luminosity. This sets our energy budget as 



AE < 



aLT 



(11) 



While we have used luminosity variations to set our energy 
budget we note that such variations will be strict upper 
limits. Indeed, it is not clear that any luminosity changes 
would be observable. If the thermal timescale of the en- 
velope is much larger than the timescale of activity cycles 
then the observed luminosity may hardly vary. To test this, 
we have examined the standard solar model (SSM) used by 
iBoothrovd fc Sackmam] (|2003D and estimated the thermal 
timescale of the convective envelope (the lower boundary of 
which we have taken to lie at a radius R = 0.713863^©). 
For each shell in the SSM the thermal energy due to fully 
ionised hydrogen, helium and associated free electrons was 
calculated (hydrogen and helium are fully ionised except for 
the very uppermost regions near the photosphere). From 
this we determine a thermal timescale of 73,000 years for 
the Sun's convective envelope. This is likely to be a lower 
limit (but a reasonable estimate nonetheless) since we have 
not included the thermal energy from metals. We conclude, 
therefore, that any luminosity variations on the timescale 
of lO's or lOO's of years may be damped considerably, and 
would be potentially unobservable. 

Armed with a prescription for defining the energy bud- 
get available to the star to drive AQ (but being aware of 
the possible reservations so far described), we can calcu- 
late orbita l perio d modulations via equation |9] Following 
lApplegatd j 19921 ) . the amplitude of the orbital period mod- 
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ulation and the amplitude of the oscillation observed in an 
O ~ C diagram are related by 
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This, combined with equation |5] and setting /3 = Ms/M, 
leads to potential O — C variations in the observed planetary 
transit times of the order, 
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where k = 2.172 x 10"''a5/3 2 . 



4 RESULTS 



lApplegcit3 l|l992f ) showed that the period changes observed 
in binaries, with just one exception, can be explained as- 
suming outer shell masses of AIs = O.IM, and that a total 
energy budget of the order of 10% of the stellar luminosity 
is available to drive the quadrupole changes. In calculating 
the orbital period changes in this work we have also adopted 
Ms = O.IM which sets P = 0.1 in equation [T31 We have 
set a = 0.1, constraining the total energy budget such that 
any RMS luminosity variations are less than 10 per cent 
of the total stellar luminosity (but note our previous com- 
ments that such luminosity variations are likely to be heav- 
ily damped). From equation 1111 this is a factor of tt lower 
than the energy budgets required by lApplegat3 lll99S ). We 
should also add that the prescription bv lLanza et al. ril998l) 
is mo re energy efficient t han that proposed by lApplegatg 
l|l992l) . The refinements bv lLanza et al.l l| 19981 ') would, there- 
fore, be able to drive larger variations than those predicted 
her e. Indeed, q u asi-pe riodic variations larger than those that 
the lApplegatg (|l992l ') prescription could sustain may have 
been observed in a few binary sys tems (see, for example, 
iBrinkworth et al]|2006l : iLanzal [20051 and references therein). 

Taking our conservative energy budget, Table[T]lists the 
expected O-C variations estimated from equation [13] for the 
parameters of 59 known transiting exoplanets. These values 
should be seen as an order of magnitude estimate of the 
Applegate effect in e xoplanet systems. St ellar rotation pe- 
riods are taken from IWatson et al.l (|2010l ) where available, 
and have been obtained from measurements of the Ca H & 
K emission lines. Other stellar rotation periods indicated by 
'~' are deduced from projected equatorial rotation veloci- 
ties, V sin i and published stellar radii from the exoplanet 
encyclopedia. Where only upper limits to vsini exist we 
can only place lower limits on the stellar rotation periods. 
Where we are unable to determine a stellar rotation period 
we have adopted 30 days in order to derive a reasonable 
estimate of the amplitude of the Applegate effect assuming 
the star has solar-like rotation. The other parameters, stel- 
lar mass, luminosity, orbital period and orbital separation 
are taken/derived from the exoplanet encyclopedia. 

The remaining variable in equation[T3]that has not been 
discussed is the modulation timescale, T, which is related to 
the period of the stellar activity cycle. For our Sun, T could 
be equal to 11 or 22 years depending on the dynamo at 
work (see Section 14.41 for a discussion). Since, in general. 



most of the host stars of exoplanets are solar-like we have 
calculated the deviations due to the Applegate mechanism 
for both of these timescales. In addition, we have also calcu- 
lated the variations over a 50 year modulation period since 
many of the orbital period variations in active binaries occur 
on longer timescales. The O — C variations derived assum- 
ing 11, 22 and 50 year modulation timescales are denoted 
5tii,5t22 and St^o, respectively. 

In addition, we have calculated the angular velocity of 
differential rotation between the outer shell and the stellar 
interior as defined in equation [4] These are expressed as the 
fractional change in angular rotation rate Af2/nii,22,50 over 
the 11, 22 or 50 year modulation period. While the addi- 
tional energy required to drive this differential rotation is 
not included in equation [TH] (since it would complicate the 
equation) we have included it in an additional code. For the 
majority of the systems in Table[l]this extra energy require- 
ment has little effect on the results. Only for those systems 
with differential rotation rates above lO"'^ does equation 1131 
start to break down. However, these still provide good 'or- 
der of magnitude' estimates for the amplitude of the O — C 
variations that may be driven by t he Applegate mec h anism 
in these systems. For comparison, lUlrich fc Bertellol (|l996[ ) 
have derived an upper limit of Sfl/Q ~ 70 x 10~^ over the 
course of the Sun's 11 year cycle. 

Inspection of Table [T] shows that the Applegate effect 
can produce large, measurable transit time variations in sev- 
eral systems. These are generally of the order of a few sec- 
onds over 11 year modulation periods, up to minutes for 
systems with 50 year stellar activity periods. Given these 
results, we now discuss the potential implications these find- 
ings have on dynamical studies of exoplanet systems. 



4.1 Implications for planet detections using TTVs 
caused by orbital perturbations by a second 
planet 

In transiting systems, the presence of a s econd planet wil l 
induce variations in the transit times ( e.g. lAgol et al.ll2005l ) 
and also durations (e.g. iKippin j|2009l ). The Applegate ef- 
fect, however, will also produce variations in the transit 
times and durations since it modulates the orbital period 
over the course of the stellar activity cycle. While the varia- 
tions due to the Applegate mechanism are at best quasi- 
periodic, there has been regular confusion as to whether 
period variations observed in some binaries are due to the 
presence o f an orbiting th ird body or the Applegate mecha- 
nism re.g. lSovdug"anll2008l ). Similarly, in exoplanet systems, 
it is possible that TTVs caused by the Applegate mecha- 
nism could mimic the perturbations caused by the presence 
of a second planet. 

To assess the possibility of confusion between TTVs 
caused by the Applegate effect and those from a bona-fide 
additional planet, we need to consider the form of the TTVs 
caused by the Applegate mechanism. First, the magnitude 
of the TTVs introduced by the Applegate mechanism drop 
off as (from equation I13p. Thus the mechanism is only 
strong for transiting planets with orbital periods of a few 
days (which encompasses most transiting systems). Second, 
the Applegate effect will cause TTVs that are modulated 
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Table 1. Summary of the estimated O-C variations induced by the Applegate mechanism for currently known transiting exoplanets. 
The first 5 columns list the system name as well as the host stars' rotation period, mass, radius and luminosity. Columns 6 and 7 list 
the planets' orbital period and orbital separation to 4 decimal places, respectively. For stars where no constraint could be placed on the 
rotation period we have assumed (for demonstrative purposes) a value of 30 days and indicated these entries with an asterisk. The next 
columns indicate the differential surface shear and expected O-C variations for stars with modulation periods of 11, 22 and 50 years, 
respectively. 



Planet 


Prot 








Porb 


a 


Af2/f2ii 






At22 


Aa/Oso 


Atso 




(days) 


(Mo) 


(R©) 


(L©) 


(days) 


(AU) 


(xlO-*) 


(s) 


(xlO-4) 


(s) 


(xlO-*) 


(s) 


WASP IQh 


10.5 


0.96 


0.94 


0.73 


0.7888 


0.0164 


1.9 


2.8 


2.7 


8.0 


4.1 


27.5 


v-'Ol U t ( U 


23.0 


0.93 


0.87 


0.53 


u.ooou 


0.0172 


1.0 


0.8 


1.4 


2.2 


2.1 


7.5 


WA Qp_1 «h 

VV Jt.Ox^ -LOU 




1.25 


1.22 


2.2* 


0.9415 




6.1 


11.5 


8.6 


32.4 


13.8 


111.0 


WA9P-1 9h 


~36.0 


1.35 


1.57 


3.49 


1.0914 


n n99Q 


1.2 


4.3 


1.7 


12.3 


2.6 


42.1 


OnT,F,-TR-'SRh 


26.3 


1.17 


1.32 


2.1* 


1.2119 


yj. KjZiiwj 


2.0 


3.0 


2.8 


8.4 


3.6 


28.7 


_L1 J_LjO O U 


^27.0 


0.924 


0.81 


0.64 


1.3062 




2.3 


0.3 


3.2 


0.9 


4.9 


3.1 


WASP dh 


^29.0 


0.9 


1.15 


1.09 


1.3382 


0.0230 


2.1 


1.6 


3.0 


4.5 


4.5 


15.3 


OCT.E-TR-1 1 


3.2 


0.78 


0.77 


0.2* 


1.4328 


0.0229 


14.9 


1.5 


21.1 


4.4 


31.8 


14.9 




'-^10.8 


0.95 


1.11 


1.39 


1.5090 


0.0254 


8.1 


3.1 


11.5 


8.9 


17.3 


30.5 


WASP-5b 


^15.7 


1.021 


1.08 


1.26 


1.6284 


0.0273 


6.1 


1.5 


8.6 


4.1 


13.0 


14.1 




30.0* 


1.26 


1.34 


2.40 


1.6899 




3.5 


1.4 


4.9 


4.0 


7.4 


13.8 




4.5 


0.97 


0.90 


0.73 


1.7430 


0.0281 


22.8 


1.9 


32.3 


5.3 


48.7 


18.2 




5.0 


1.24 


1.31 


2.59 


1.8468 


0.0317 


26.8 


7.8 


38.0 


22.1 


57.2 


75.9 


WA 9P-9h 


30.0* 


0.84 


0.83 


0.46 


2.1522 


0.0314 


4.8 


0.2 


6.8 


0.5 


10.3 


1.6 


HAT- P- 7b 


^25.0 


1.47 


1.84 


4.95 


2.2047 


0.0379 


6.8 


4.5 


9.7 


12.8 


14.6 


43.8 


nLj ±o,j 1 oou 


13.2 


0.80 


0.79 


0.34 


2.2186 


0.0310 


12.0 


0.3 


16.9 


0.8 


25.5 


2.7 


WASP-14h 


'-^13.4 


1.319 


1.30 


2.66 


2.2438 


0.0370 


9.8 


1.9 


13.8 


5.3 


20.9 


18.3 


TrF,S-2h 


24.8 


0.98 


1.00 


1.05 


2.4706 




9.0 


0.4 


12.7 


1.1 


19.2 


3.7 


OnT,F2-TR-T,Qh 


~2.0 


1.52 


1.53 


4.86 




o'^ns 

V J . VJ ij V./ C.J 


129.4 


38.1 


183.1 


107.8 


276.0 


369.4 


WASP-lb 


30.0* 


1.24 


1.38 


2.54 


2.5200 


0382 


7.7 


1.1 


10.9 


3.1 


16.5 


10.6 


XO-2b 


>21.0 


0.98 


0.96 


0.68 


2.6158 


n (1369 


9.9 


<0.3 


14.1 


<0.8 


21.2 


<2.8 




48.0 


0.452 


0.46 


0.04 


2.6439 


0.0287 


3.3 


0.0 


4.6 


0.0 


7.0 


0.1 


HAT-P-'ib 


~22.7 


1.16 


1.17 


1.55 


2.7885 


0.0408 


12.0 


0.6 


16.9 


1.6 


25.5 


5.4 




~12.6 


1 ^ 


1.50 


2.88 


2.8759 


0.0431 


23.1 


2.8 


32.6 


7.9 


49.1 


27.0 


HAT-P-Sb 


~80.0 


0.936 


0.82 


0.44 


2.8997 




3.1 


0.0 


4.4 


0.1 


6.6 


0.3 


HAT-P-1 "Vn 


~27.0 


1.22 


1.56 


2.21 


2.9163 


0.0426 


9.6 


1.5 


13.6 


4.3 


20.5 


14.7 


TrES- lb 


33.5 


0.87 


0.82 


0.46 


3.0301 


VJ . uoyo 


8.6 


0.1 


12.1 


0.2 


18.3 


0.8 


HAT-P-4b 


~14.6 


1.26 


1.59 


2.68 


o . uouo 


VJ . U'rt'iVJ 


20.7 


2.9 


29.3 


8.2 


44.2 


28.0 


H AT-P-Sh 


^^6.2 


1.28 


1.58 


3.32 


3.0763 


0.0487 


55.0 


6.0 


77.7 


17.1 


117.2 


58.6 


WARP-lOb 


>6.6 


0.71 


0.78 


0.26 




0.0371 


39.6 


<0.4 


55.9 


<'! 1 


84.3 


<3.9 


onT.F^TR-inh 


15.8 


1.18 


1.16 


1.32 


3.1013 


0.0416 


19.6 


0.7 


27.7 


1.9 


41.7 


6.5 


WARP-lfib 


~16.0 


1.022 


0.95 


0.76 


3.1186 


0.0421 


19.6 


0.3 


27.7 


0.8 


41.7 


2.6 


XO-3b 


~3.8 


1.213 


1.38 


2.91 


3.1915 


0.0454 


106.5 


6.6 


150.6 


18.8 


227.0 


64.4 


HAT-P-1 9h 


30.0* 


0.73 


0.70 


0.21 


3.2131 


0.0384 


9.3 


0.0 


13.2 


0.1 


19.9 


0.4 


WASP-fih 


'-^31.5 


0.88 


0.87 


0.60 


3.3610 


0.0421 


12.0 


0.1 


17.0 


0.3 


25.6 


1.1 




14.9 


1.01 


1.15 


1.47 


3.5247 


0.0471 


31.0 


0.7 


43.8 


2.0 


66.0 


6.9 




~10.8 


1.384 


1.81 


4.08 


o . oooy 


VJ . uouy 


39.2 


5.4 


55.4 


15.3 


83.5 


52.3 


Oni.F TR 211 h 


30.0* 


1.33 


1.64 


3.87 


3.6772 


0.0510 


16.6 


1.3 


23.4 


3.8 


35.3 


13.1 


WA SP lib 


~27.3 


0.82 


0.81 


0.36 


3.7225 


0.0439 


14.5 


0.1 


20.5 


0.2 


30.9 


0.8 


WASP-1 7b 


~7.8 


1.2 


1.38 


3.15 


3.7354 


0.0510 


74.2 


2.7 


104.9 


7.7 


158.2 


26.5 


WA9P-1 


~18.7 


1.18 


1.48 


3.09 


3.7521 


0.0542 


29.1 


1.3 


41.2 


3.8 


62.1 


13.1 


WAT' P f\K 




1 OQ 


1 Afi 


O.O 1 


o.oOoU 




1 U.O 


9 Q 


QQ 4 




J-OU.U 


97 R 


Lupus-TR-3b 


30.0* 


0.87 


0.82 


0.38 


3.9141 


0.0464 


14.5 


0.1 


20.6 


0.2 


31.0 


0.6 


HAT-P-9b 


~5.5 


1.28 


1.32 


2.55 


3.9229 


0.0530 


105.7 


2.5 


149.2 


7.0 


225.3 


23.8 


XO-lb 


>16.0 


1.00 


0.93 


1.0* 


3.9415 


0.0488 


37.0 


<0.2 


52.3 


<0.6 


78.8 


<2.1 


OGLE-TR-182b 


30.0* 


1.14 


1.14 


1.44 


3.9791 


0.0510 


18.4 


0.2 


26.0 


0.7 


39.2 


2.3 


OGLE-TR-lllb 


30.0* 


0.82 


0.83 


0.41 


4.0145 


0.0470 


16.1 


0.1 


22.8 


0.2 


34.4 


0.7 


Corot-5b 


~90.0 


1.00 


1.19 


1.75 


4.0379 


0.0495 


7.1 


0.1 


10.1 


0.4 


15.2 


1.3 


XO-4b 


~8.9 


1.32 


1.55 


2.28 


4.1250 


0.0555 


57.3 


2.4 


81.0 


6.7 


122.1 


23.0 


XO-5b 


~80.0 


0.88 


1.06 


0.93 


4.1878 


0.0487 


7.5 


0.1 


10.6 


0.2 


16.0 


0.9 


Corot-3b 


~4.6 


1.37 


1.56 


4.52 


4.2568 


0.0570 


162 


6.0 


229.1 


16.9 


345.4 


57.8 


WASP-13b 


>13.8 


1.03 


1.34 


1.86 


4.3530 


0.0527 


48.6 


<1.2 


68.8 


<3.5 


103.7 


<12.1 


HAT-P-lb 


19.7 


1.133 


1.12 


1.42 


4.4653 


0.0553 


35.9 


0.3 


50.7 


0.9 


76.5 


3.0 


HAT-P-llb 


30.0* 


0.81 


0.75 


0.26 


4.8878 


0.0530 


21.2 


0.0 


30.0 


0.1 


45.3 


0.3 


WASP-7b 


~3.7 


1.28 


1.24 


2.30 


4.9547 


0.0618 


2.54.1 


2.0 


359.4 


5.5 


541.8 


19.0 


HAT-P-2b 


4.0 


1.36 


1.64 


3.79 


5.6335 


0.0689 


282.1 


5.3 


398.9 


14.9 


601.4 


51.0 


Corot-6b 


30.0* 


1.055 


1.03 


1.30 


8.8866 


0.0855 


100.7 


0.0 


142.4 


0.2 


214.7 


0.6 


Corot-4b 


8.9 


1.1 


1.15 


1.75 


9.2021 


0.0900 


369.8 


0.3 


523.0 


0.9 


788.5 


3.1 


HD 17156b 


22.1 


1.24 


1.45 


2.57 


21.2169 


0.1623 


715.0 


0.1 


1011.2 


0.3 


1524 


1.0 



6 C.A. Watson, and T. R. Marsh 



over the timescale of the steUar activity cycle, most likely 
over the course of years to decades. The question then be- 
comes, are there any possible mechanisms by which a second 
planet can modulate the transit times of short period planets 
by a few seconds over the course of ~10 years? 

lAgol et al.1 (|2005h outline a number of mechanisms by 
which a second planet can cause TTVs. Confusion with inte- 
rior perturbing planets can be immediately rejected since the 
modulation of the orbital period will occur on far too short 
a timescale, in order to cause TTVs modulated on decade 
timescales requires an exterior planet on a relatively wide 
orbit. Furthermore, TTVs caused by an additional planet 
in a mean-motion resonance orbit generally drive signifi- 
cantly larger transit-time variations than those possible by 
the Applegate mechanism. However, for exterior planets on 
eccentric orbits with much l arger orbita l periods (and not 
in mean-motion resonance), lAgol et"al] |2005) found that 
deviations in the transit times accumulate over the orbital 
period. Pout, of the outer planet to give, 

St = iiouteout \ J Pout (14) 

\ flout / 

where flout is the star/planet mass ratio for the outer planet, 
Gout is the eccentricity of the outer planet, and Uin and aout 
are the orbital separations of the inner and outer planets, 
respectively. 

Figure [1] shows the magnitude, 5t, of the transit tim- 
ing variations caused by the Applegate mechanism (solid 
lines) compared to gravitational perturbation by an exte- 
rior planet (dashed line) as a function of the orbital radius, 
ain, of the interior transiting planet. For the Applegate de- 
viations we have taken the example of a solar-like host star 
(M* = IMq, R, = li?0, L» = li©) and, as before, set the 
parameters a — 0.1 and /3 = 0.1 and assumed a modulation 
timescale, Pmod, of 11 years. We have then plotted three 
curves corresponding to stellar rotation rates of 5, 10 and 
20 days. To calculate the maximum deviations expected as 
a result of gravitational perturbations by a second planet 
we have assumed a 13A/j exterior planet (the accepted up- 
per mass-limit for a planet) in a highly eccentric [eout = 
0.9) orbit (dashed line in Fig. [T]) with an orbital period of 
11 years (matching the modulation period assumed in the 
Applegate calculation). This, therefore, represents the max- 
imum deviations expected from a second planet. For com- 
parison, a Jovian-mass object in a Jupiter-like orbit would 
cause deviations with an amplitude l/208th that of the ob- 
ject considered here. 

Fig.[T]shows that, for even modest stellar rotation rates, 
the Applegate mechanism is likely to dominate over TTVs of 
short period (a < 0.03AU) transiting systems around solar- 
like stars caused by the gravitational perturbations of a sec- 
ond planet on a long period, non-resonant orbit. For modula- 
tion timescales longer than 11 years (corresponding to longer 
activity cycles or longer exterior planet orbital periods) the 
crossover points of the curves occurs at larger inner-planet 
orbital separations, scaling as ain^xover oc P^lod- Since the 
Applegate mechanism becomes less efficient as the stellar 
mass (and luminosity) decreases, while gravitational per- 
turbations due to a second planet will increase (due to the 
flout factor in equation 1 14|) this situation is likely to reverse 
for planetary systems around lower-mass host stars. Thus, 
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Figure 1. The magnitude, &t, of the transit timing variations 
caused by the Applegate mechanism (solid lines) compared to 
gravitational perturbation by an exterior planet (dashed line) as 
a function of the orbital radius, ain, oi the interior transiting 
planet. The deviations caused by the Applegate mechanism as- 
sume a Sun-like host star with rotation rates (indicated on the 
plot) corresponding to a 5, 10 and 20 day stellar rotation period, 
respectively, and a modulation period Pmod = H yrs. The dashed 
line represents the TTVs caused by a 13Mj mass planet with an 
11 year orbital period [aout = 4.95AU) on a highly eccentric {eout 
= 0.9) orbit. 

searches for TTVs around M-dwarf stars may not only be 
more sensitive to the presence of planetary companions, but 
may also be less prone to 'confusion' with TTVs caused by 
the Applegate effect. (In addition, M-dwarfs may not have 
an appreciable dynamo mechanism to drive the Applegate 
effect). 

We have not considered TTVs for planets in mean- 
motion resonance orbits which, under the right conditions, 
can drive significantly larger deviations than those described 
above. The conditions for possible confusion with the Ap- 
plegate mechanism (we require long modulation timescales 
of years to decades and a short orbital period for the tran- 
siting planet) restricts the perturbing exterior planet to long 
orbital periods. The large difference between the orbital pe- 
riods of the interior and exterior planets means that there 
will be no strong resonance between the components. Mean 
motion resonances are unlikely to provide significantly in- 
creased gravitational perturbations si nce the examples co n- 
sidered here have period ratios >150 (|Beauge et al.|[2005h . 

4.2 Implications for planet detections using TTVs 
caused by light travel time effects 

The presence of a second planet in a long period may also 
cause TTVs through the variation of the light travel time 
due to the refiex motion of the host-star. Indeed. iDeeg et ahl 
l|2000D claim tentative evidence for a Jovian mass object in 
the eclipsing binary CM Dra by observing this effect. For an 
exterior planet of mass Mp on an orbit with a semi-major 
axis aout, the amplitude of the timing deviation is. 
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for a solar-mass star. These periodic variations naturally 
occur over a timescale equal to the orbital period of the 
outer planet. 

For a Jovian-mass object in a Jupiter-like orbit this 
amounts to a periodic deviation of ~5 seconds over the 
course of ~12 years. This is of a similar magnitude and 
timescale as those potentially driven by the Applegate ef- 
fect. There is, therefore, considerable scope for confusion 
between TTVs caused by light travel time delays due to the 
presence of a second planet and those caused by the Apple- 
gate effect. The only secure means of distinguishing between 
these two effect is to look for strictly periodic variations since 
the Applegate mechanism will be quasi-periodic at best. 



4.3 Implications for measuring tidal dissipation 

Recently a subset of exoplanets orbiting extremely close to 
their host stars with periods less than 1 day have been dis- 
covered. Prime examples of these are the planets WASP- 18b, 
whic h was dis covered to be orbiting a 1.24M0 star every 0.94 
days (iHellier et a l., 2009 ): WASP-19b (orbital period = 0.79 



Hebb et al 



Leger et al 



and Corot-7b (orbital period : 



0.85 



20091 ). These planets are subject to large 



days; 
days; 

tidal forces and can be used as tests of tidal dissipation the- 
ory. In this scenario, the close proximity of the exoplanet 
to the host star raises a tidal bulge on the stellar surface. 
This in turn exerts an additional torque which drains angu- 
lar momentum from the planetary orbit for systems where 
the orbital period is shorter than the stellar rotation period. 

In the case of WASP-18, the tidal dissipation is such 
that assuming a tidal quality factor D (this is normally de- 
noted Q but we have changed this to distinguish it from the 
quadrupole moment used earlier) found from studies of bi- 
naries and the giant planets in our Solar System means that 
the planet would be tida lly disrupted i n ^1 Myr. Since the 
system's age is ~5 Gyr, iHellier et al.1 (|2009l 'l conclude that 
either they have caught WASP-18b in an e xtremely rare 
state, or that D is much higher than expected. IHellier et al.l 
l|2009h quote a period change of -0.00073 (10® /D) s yr'^ 
For our solar system, D = 10^ — 10® leading to orbital pe- 
riod changes of 28 s after 10 yr for D = 10®. Therefore, 
transit timing studies of WASP- 18b are vitally important 
since the tidal decay should be directly measurable in this 
system. This would allow our understanding of tidal dissi- 
pation to be tested and would also allow the stellar interior 
to be probed since the value of D depends on how waves 
generated by tides are dissipated. 

The magnitude of the Applegate effect we have esti- 
mated for WASP-18 leads to 0-C variations of the order 
of lO's of seconds on decade timescales - approximately the 
same as predicted due to tidal decay. Such variations could 
easily either mask, or mimic (briefly), or even temporarily 
reverse any period change caused by tidal decay, especially if 
the tidal quality factor is higher than 10®. Thus, for WASP- 
18b, caution has to be exercised in interpreting the nature of 
any detected transit-time variations, since these will consist 
of a quasi-periodic variation on the timescale of the stellar 
activity cycle super-imposed on a steady trend of decreas- 
ing period due to orbital decay. Only by observing over the 
course of an activity cycle could one begin to reliably dis- 
tinguish between the two mechanisms. Thus, WASP-18 is 



unlikely to reveal its tidal dissipation history (or future) for 
several years yet. 

4.4 Implications for stellar dynamo theory 

Finally, we discuss the use of accurate transit timing as a 
means to probe the nature of the magnetic field generat- 
ing dynamo operating in the host stars' interior. The study 
of the Applegate mechanism via eclipse timings of eclips- 
ing binaries is often confused with orbital period variations 
caused by angular momentum losses or exchanges. These 
include effects due to magnetic braking and mass transfer. 
We propose that transiting exoplanet systems, by compar- 
ison, provide much cleaner laboratories in which to study 
such effects. Furthermore, unlike in the tidally locked binary 
systems, exoplanet systems have the crucial ingredient of a 
difference between the stellar rotation period and its orbital 
period (the orbital period of the exoplanet). If variations 
are seen at a level that correlates with the predictions of 
Table [T] for two stars with similar orbital periods but differ- 
ent rotation periods, then this would provide good evidence 
that the Applegate mechanism is at work. Naturally, how- 
ever, one should be aware that gravitational perturbations 
caused by an additional long-period planet in the system 
could initially be confused with variations due to the Ap- 
plegate mechanism. Thus, the quasi-periodic nature of such 
variations has to be firmly established before one can be 
confident that such TTVs can be used to probe the stellar 

dyn amo. 

iLanza et al.l (|l998l ) showed that different dynamo mech- 
anisms result in different observable manifestations of the 
Applegate effect. If an ail dynamo is in operation then the 
quasi-cyclic orbital period modulations should occur on the 
same timescale as the spot coverage modulation (e.g. over 
11 years in the case of our Sun). In contrast, if an a^Q 
dynamo operates then the observed orbital period modu- 
lation should occur over a timescale twice as long as the 
period of the spot modulation (~ 22 years in the case of 
our Sun). Given the large number of transiting exoplanet 
systems with host stars (with widely varying fundamental 
parameters such as age, rotation rate, masses etc.) that will 
be discovered by space-missions such as Kepler and Plato, 
long-term systematic monitoring for TTVs may reveal the 
nature of the dynamo mechanism at work. 



5 DISCUSSION 

We have estimated the likely transit timing variations in- 
duced by changes in the quadrupole moment of the host 
star in transiting exoplanet systems driven by the Apple- 
gate effect. Depending on the length of the activity cycle, 
TTVs of several minutes are plausible for a number of the 
currently known transiting exoplanets. While the timescales 
and sizes of Applegate driven TTVs are of the wrong magni- 
tude to be confused with TTVs driven by additional planets 
in mean motion resonances, there appears to be much scope 
for confusion with TTVs caused by light travel time effects 
caused by massive, Jupiter-like planets on wide orbits. 

The magnitude of TTVs driven by the Applegate mech- 
anism also grow as the star-planet separation decreases (as- 
suming all other factors are equal). Indeed, for the shortest 
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period transiting exoplanets, such as WASP-18b, the Ap- 
plegate mechanism could potentially be mistaken for orbital 
period changes due to tidal dissipation. Indeed, the orbital 
decay due to tidal dissipation could even be temporarily re- 
versed since the orbital period variations due to the Apple- 
gate effect can take on either sign. We therefore urge caution 
when interpreting TTVs, especially those that appear to be 
occurring on timescales of years to decades. In all cases, the 
clear signature that the Applegate effect is at work is that 
the TTVs are quasi-periodic. Only once the strict periodic- 
ity of any TTV has been ascertained should investigators be 
confident in their final interpretation. 

Finally, an alternative mechanism for driving quasi- 
periodic stellar quadrupole variati ons (and hence orbi tal pe- 
riod variations) was put forward bv lLanza etall l|l998h . This 
rests on the principle that changes in the azimuthal mag- 
netic field intensity can change the stellar quadrupole mo- 
ment by altering the effective centrifugal acceleration, the 
stellar dynamo effectively interchanges magnetic energy and 
rotational kinetic energy. The main feature of note in this 
prescription is that relative changes in the angular veloc- 
ity required to drive orbital peri od changes are a factor of 
2 smaller than that required by lApplegeit3 (|l992l ). While 
it i s difficult to assess what TTVs may be expected from 
the iLanza et al.l |l998) work, it is quite possible that non- 
periodic TTVs with magnitudes exceeding those outlined in 
this paper could be observed. This latter point should also 
be taken in light of the fac t that we have been very con- 
servative (in comparison to lApplegate|[T993 ) in the energy 
budget we have assumed to be available for driving the stel- 
lar quadrupole variations. Indeed, as outlined ear lier, quasi- 
period ic variations larger than those that the lApplegatg 
l| 19921 ) prescription could sustain may already have been 
observed in a few binary systems. Studies of the long-term 
trends in the transit times of short-period exoplanets could 
provide crucial evidence for settling may of the points above, 
presenting valuable insights into the working of stellar dy- 
namos across a broad range of fundamental stellar proper- 
ties. 



ACKNOWLEDGMENTS 

This research has made use of NASA's Astrophysics Data 
System Bibliographic Services. The authors also gratefully 
acknowledge the use of the Extrasolar Planets Encyclopedia 
|http://exoplanet.eu/[ ). 



REFERENCES 

Agol E., Steffen J., Re'em S., Clarkson W., 2005, MNRAS, 
359, 567 

Applegate J. H., 1992, ApJ, 385, 621 

Applegate J. H., Patterson J., 1987, ApJ, 322, L99 

Baptista R., Borges B. W., Bond H. E., Jablonski P., 

Steiner J. E., Grauer A. D., 2003, MNRAS, 345, 889 
Baptista R., Catalan M. S., Costa L., 2000, MNRAS, 316, 

529 

Baptista R., Jablonski F., Oliveira E., Vrielmann S., Woudt 
P. A., Catalan M. S., 2002, MNRAS, 335, 75 



Baptista R., Jablonski F., Steiner J. E., 1992, A&A, 104, 
1557 

Bean J. L., 2009, A&A, 506, 369 

Beauge C, Callegari J., Ferraz-Mello S., Michtchenko 
T. A., 2005, in Knezevic Z., Milani A., eds, lAU CoUoq. 
197: Dynamics of Populations of Planetary Systems, Pro- 
ceedings Resonances and stability of extra-solar planetary 
systems. Cambridge University Press, Cambridge, p. 3 
Bond I. A., Freeth R. V., 1988, MNRAS, 232, 753 
Boothroyd A. I., Sackmann I.-J., 2003, ApJ, 583, 1004 
Borges B. W., Baptista R., Papadimitriou C, Giannakis 

O., 2008, A&A, 480, 481 
Brinkworth C. S., Marsh T. R., DhiUon V. S., Knigge C, 

2006, MNRAS, 365, 287 
Deeg H. J., Doyle L. R., Kozhevnikov V. P., Blue J. E., 

Martin E. L., Schneider J., 2000, A&A, 358, L5 
Echevarn'a J., Alvarez M., 1993, A&A, 275, 187 
Gibson N. P., PoUacco D., Simpson E. K., Barros S., Joshi 
Y. C, Todd I., Keenan F. P., SkiUen I., Benn C, Christian 
D., Hrudkova M., Steele I. A., 2009, ApJ, 700, 1078 
Gibson N. P., PoUacco D. L., Barros S., Benn C, Christian 

D. , Hrudkova M., Joshi Y. C, Keenan F. P., Simpson 

E. K., SkiUen I., Steele I. A., Todd I., 2010, MNRAS, 401, 
1917 

Glownia Z., 1986, Adv. Space Res., 6, 211 

HaU D. S., Kreiner J. M., 1980, Acta Astronomica, 30, 387 

Hebb L., Collier-Cameron A., Triaud A. H. M. J., Lister 
T. A., Smalley B., Maxted P. F. L., Hellier C, Anderson 
D. R., PoUacco D., Gillon M., Queloz D., West R. G., 
Bentley S., Enoch B., Haswell C. A., Home K., Mayor 
M., Pepe F., et al. 2010, ApJ, 708, 224 

Hellier C, Anderson D. R., Cameron A. C, Gillon M., 
Hebb L., Maxted P. F. L., Queloz D., Smalley B., Triaud 
A. H. M. J., West R. G., Wilson D. M., Bentley S. J., 
Enoch B., Horne K., Irwin J., Lister T. A., Mayor M., 
Parley N., Pepe F., et. al 2009, Nature, 460, 1098 

Kipping D. M., 2009, MNRAS, 392, 181 

Lanza A. F., 2005, MNRAS, 364, 238 

Lanza A. F., Rodono M., Rosner R., 1998, MNRAS, 296, 
893 

Leger A., Rouan D., Schneider J., Barge P., Fridlund M., 
Samuel B., OUivier M., Guenther E., DeleuU M., Deeg 
H. J., Auvergne M., Alonso R., Aigrain S., Alapini A., 
Almenara J. M., Baglin A., Barbieri M., et al. 2009, A&A, 
506, 287 

Marsh T. R., Pringle J. E., 1990, ApJ, 365, 677 
Mayor M., Queloz D., 1995, Nature, 378, 355 
Miller-Ricci E., Rowe J. F., Sasselov D., Matthews J. M., 
Guenther D. B., Kuschnig R., Moffat A. F. J., Rucinski 
S. M., Walker G. A. H., Weiss W. W., 2008, ApJ, 682, 
586 

Miralda-Escude J., 2002, ApJ, 564, 1019 

Pal A., Kocsis B., 2008, MNRAS, 389, 191 

Rabus M., Deeg H. J., Alonso R., Belmonte J. A., Almenara 

J. M., 2009, A&A, 508, 1011 
Rafikov R. R., 2009, ApJ, 700, 965 
Saar S. H., Brandenburg A., 1999, ApJ, 524, 295 
Simon A., Szatmary K., Szabo G. M., 2007, A&A, 470, 727 
Soydugan F., 2008, Astron. Nachr., 329, 587 
Ulrich R. K., Bertello L., 1996, ApJ, 465, L65 
Warner B., 1988, Nature, 336, 129 

Watson C. A., Littlefair S. P., Collier Cameron A., DhiUon 



Orbital period variations of hot-Jupiters 9 

V. S., Simpson E. K, 2010, MNRAS, submitted 
Wolf S., Mantel K. H., Home K., Barwig H., Schoembs R., 
Baerbantner O., 1993, A&A, 273, 160 



